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fireball. If the inverse Compton process dominates the cooling of shocked elec- 
trons, an X-ray flare produced by the first order Compton scattering would 
emerge in the very early afterglow phase, with the bulk of the shock energy 
radiated in the second order scattering component at 10-100 MeV. The domi- 
nance of inverse Compton cooling leads to the lack of prompt optical flashes. We 
show that for plausible parameters this scattering process can produce an X-ray 
flare with a relative amplitude change by a factor of several. Flares with a larger 
; amplitude and multiple X-ray flares in a single event are likely to be produced 

by another mechanism (e.g. internal shocks). 
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1. Introduction 



The Swift satellite is a multi-wavelength observatory designed to detect GRBs and 
^ ■ their X-ray and UV/optical afterglows. Thanks to its fast pointing capabilities Swift is 
disclosing the early afterglow phase. The Swift X-Ray Telescope (XRT) found that most 
X-ray afterglows fall off rapidly for the first few hundred seconds, followed by a less rapid 
decline (Tagliaferri et al. 2005). In the early afterglows of GRB 050406 and GRB 050502b, 
XRT detected mysterious strong X-ray flares: rapid brightening of the X-ray afterglow after 
a few hundred seconds post-burst (Burrows et al. 2005). These results suggest the existence 
of additional emission components in the early afterglow phase besides the conventional 
forward shock emission. 
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In recent years, ground-based robotic telescopes reported the lack of prompt optical 
emission, except for a few cases. The Swift UV/Optical Telescope (UVOT) provided further 
stringent upper limits at very early epochs after the bursts (Roming et al. 2006). The limits 
~ 20 mag at < 100 sec are much lower than the 9th magnitude optical flash associated with 
GRB 990123 (Akerlof et al. 1999). Although the host extinction at z ~ 1 and/or the high 
magnetization of a fireball (Zhang & Kobayashi 2005) can explain the lack of optical flashes, 
other suppression mechanisms might be involved. 

According to the standard relativistic fireball model, reverse shocks are expected to 
radiate photons in the optical/IR band via the synchrotron process in a very early afterglow 
phase (Meszaros & Rees 1997; Sari & Piran 1999). Although the contribution of the reverse 
shock synchrotron emission to the X-ray band is small, electrons in the reverse-shocked 
region can up-scatter the synchrotron photons (synchrotron self-inverse Compton emission, 
hereafter SSC) to the X-ray or even higher energies 1 (Meszaros & Rees 1993; Wang, Dai & 
Lu 2001a,b; Granot & Guetta 2003). The SSC emission should produce additional features 
in the early X-ray afterglows. If inverse-Compton (IC) cooling is the dominant cooling 
mechanism of the electrons, the energy available for the synchrotron process is significantly 
reduced, and prompt optical flashes could be faint (Beloborodov 2005). In this paper, we 
study the role of the reverse shock SSC emission in early X-ray afterglows. 



2. Reverse Shocked Ejecta 

We consider a relativistic shell (fireball ejecta) with an isotropic energy E and an initial 
Lorentz factor T expanding into a homogeneous ambient medium with particle number 
density n. The evolution of reverse shocks is classified into two cases (Sari & Piran 1995) by 
using a critical Lorentz factor T c = (3(1 + zfE /327rnm p c 5 T 3 ) 1/8 - 130 C 378 ^ 8 ^' 378 ^ 178 
where C = (1 + z)/2, E 52 = E/10 52 ergs, T 2 = T/10 2 sec and T is the duration of the GRB. 
If T > T c (thick shell case), the shell is significantly decelerated by a reverse shock. The 
Lorentz factor at the shock crossing time t d ~ T is given by T d ~ T c . If T < T c (thin shell 
case), the reverse shock can not decelerate the shell effectively and ~ To- The deceleration 
time t d = (l + z)(3E/327vTlnm p c 5 ) 1/3 is larger than the GRB duration. 

We assume that constant fractions e e and e# of the shock energy go into the electrons and 



1 Previous studies mainly address the prompt high-energy gamma-rays (> 100 MeV) detected by EGRET. 
Higher order scattering was ignored. As we will show, for plausible parameters the second order IC is 
important to discuss early afterglow. The correction significantly modifies the amplitude of an X-ray flare 
(the first order IC component) when the Compton parameter is large. 
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magnetic fields, respectively, and that electrons are accelerated in the shock to a power-law 
distribution Nd^ oc 7^7,7 > lm ~ [(p-2)/(p-l)]e e (m p /m e )(T /T d ) ~ 180(e e /0.3)(r /r d ) 
where p = 2.5 was assumed. Using the number of electrons in the shell N e and the decel- 
eration radius Rd ~ 2cr^/(l + 2), the optical depth of the shell is r(Rd) = <JTN e / AnR 2 d = 
(l/3)lZMO'TnRd where we used the fact that the mass of the shell is larger by a factor of 
T^m = r^/r than that of the ambient material swept by the shell at the deceleration time. 



3. The First and Second order IC 

Reverse-shocked electrons emit optical photons via the synchrotron process. Since the 
random Lorentz factors of electrons are typically a few hundreds, the up-scattered photons 
are in X-ray band. If the first order IC is in the X-ray band (hux ~ 5 keV), the comoving 
photon energy of the first IC is hv' x = hvx/F ~ 50 eV where T is the bulk Lorentz factor of 
the shocked shell. The Klein-Nishina effect does not suppress the second order IC scattering 
as long as the random Lorentz factor 7 m is below ~ m e c 2 /hv' x ~ 100r. The second IC 
component appears in 10 — 100 MeV range. Higher order scattering (three or more) can be 
ignored because of the Klein-Nishina effect. 

The ratio of the IC to synchrotron luminosity can be computed in a general way (Sari 
& Esin 2001). The luminosity ratios, in the limit of up to second scattering, are given by 
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where U syn , Uic\i s t, Uic,2nd, Ub and U e are the energy density of synchrotron radiation, the 
1st IC, 2nd IC, magnetic field and random electrons, respectively. We used r]U e = U syn + 
Uic,ist + Uic,2nd- V is the fraction of the electron energy that was radiated away: r\ = 1 for 
fast cooling and r\ = (7 c /7m) 2 ~ p f° r sl° w cooling. Solving eq (1) for x we obtain 

_ f (W e .B) if (W^b) < 1 ( as 

\ ( V e e /e B y/ 3 if (rje e /6 B ) » 1 {6) 

When the IC emission dominates the overall electron cooling, it reduces the energy available 
for synchrotron radiation. Consequently, the cooling-break energy of the electron distribution 
7 C is reduced from the synchrotron-only value 7 CjS by a factor of (1 + x + x 2 ). 

If (e e /e B ) 3> 1 and r] s = ( / y c ,s/lm) 2 ~ p <sC 1, the efficiency rj depends on x. From eq. (1), 
one obtains x ~ (i] s e e /e B ) for (e e /e B ) < r\, x ~ (^ s e e / 'cb) 1 ^ 7-2 ^ for r\ < (e e /e B ) < r 2 , and 



-4- 



x ~ (ee/es) 1 ^ 3 for (e e /eB) > r 2 where ri = rjj 1 and r 2 = rj s 3 . A smaller e# makes 
the Compton parameter x larger, but the dependence is rather weak. Considering scalings 
i] s oc e p B 2 , ri oc t B p ~ e^ 2 and r 2 oc e B 3 ^ 2 , the second regime r x < (e e /es) < r 2 should be 
achieved in the limit of e# — > 0. Then, we obtain x oc e^ 3 p ^ 7 2p ^ ~ e B^ 4 for P = 2.5. 



4. SSC Emission from Reverse Shock 

The synchrotron spectrum of reverse shock emission can be approximated by a broken 
power-law with break frequencies, typical frequency v m and cooling frequency v c (e.g. Sari, 
Piran & Narayan 1998). The 1st IC spectrum is also roughly described by a broken power 
law with break frequencies v 1 ^ and v^F . The spectral characteristics of the 1st SSC emission 
are given by 

v m ~ 2 7 ^z/ m , vf ~ 2 7c 2 z/ c , ~ /trF ma;E , (4) 

where k is a correction coefficient. The ratio of vF v peaks also gives their luminosity ratio. 
It is L IC ,ist/L syn ~ v>Fl c (VpC ak )/vF u (i>p eak ) = 2KT^ m ^ c r] where v peak = max[z/ m , v c ] and 
Vpeak = max [ 1/ m j v c C ]- We can reduce this estimate to eq. (3) with a normalization k = 
4(p-l)/(p-2). 

To assess the relative importance of the 1st IC emission and the synchrotron emission 
in the X-ray band (the 2nd IC component is well above X-ray band), we consider the flux 
ratio at the IC (F v ) peak. For slow cooling, the IC spectrum peaks at i/f , and declines as 
Fl° oc //-(p^ 1 )/ 2 for z/f < v < z/ C an d as ^1° v~ p ^ 2 above z/ c , while the synchrotron 
spectrum peaks at v m , and declines as F u oc z/-^ 1 )/ 2 for v m < v < v c and as F v oc z/~ p / 2 
above z/ c . The flux ratio at the Fj c peak is given by 

F IC( ic\ F ic 

v Y_jn J _ * max < ~p-2 /r\ 

where p < 3 was assumed. The effective Compton parameter x = 2/tr7 rn 7 c r]; the equality is 
achieved only when v c ~ z/ m . If electrons are in the fast cooling regime, the 1st IC peaks at 
h>l c . We can show that the flux ratio at v* c has the same upper limit < ^^x. Therefore, 
the contrast of a SSC bump should be less than ^ m ~ 2 x ~ lOx (e e /0.3) p ~ 2 (r /T rf ) p ~ 2 where 
(r /r d ) ~ 1 for a thin shell. The bump could be more significant in the thick shell case, but 
in this case the SSC emission occurs at the end of the prompt gamma-ray emission, and it 
is difficult to separate the SSC emission from the internal shock signal. 

If the 6b parameter is small, the IC cooling becomes more important compared to the 
synchrotron cooling. However, as we showed at the end of the previous section, the depen- 
dence is weak: x oc e B ^ 4 . To achieve x >> 1, a very small e# is required. The cooling 
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frequency v c then becomes much higher than the typical frequency v, m . It is possible to show 
that the flux ratio at z/f is insensitive to e B as F/ c / 'F v {y^) = {F^J F max )(u^ /v m )^l 2 ~ 
ary^ -3 oc e^. For plausible parameters, x could be several at most. If x is large, most energy 
is radiated in the 2nd IC component, the SSC bump and the base line X-ray afterglow should 
be faint. The synchrotron emission (optical flash) should be highly suppressed. 

5. Forward Shock Emission 

A forward shock also emits X-rays via the synchrotron process. At the deceleration 
time, the spectral characteristics of the forward and reverse shock synchrotron emission are 
related (e.g. Zhang, Kobayashi & Meszaros 2003) as 

v m ,f ~ Tl l ^ 2 "R 2 M v rn) v c j ~ TZ B 3/2 1Z x 2 u c , F max j ~ 'R B l2 ll~^ F max (6) 

where the subscript / denotes forward shock quantities, we have assumed that the e e and p 
parameters are the same for both shock regions, but with different e# as parameterized by 
TZg = (cbj/cb)- The reason we introduce the 1Z& parameter is that the fireball ejecta could 
be endowed with a primordial magnetic field (e.g. Zhang et al. 2003 and references therein). 
TZx = (1 + Xf)/(1 + x + x 2 ) is a correction factor for the IC cooling. In the forward shock 
region, a once-scattered synchrotron photon generally has energy larger than the electron 
mass in the rest frame of the second scattering electrons. Multiple scattering of synchrotron 
photons can be ignored. We obtain Xf = (i]e e /eBj) for (r)e e /eBj) <§C 1 and Xf = {r\e e j 'e^,/) 1 / 2 
for (r)e e /eB,f) S> 1 (Sari & Esin 2001). TZm = r^/r is the mass ratio defined in section 2. 

The 1st IC peak is much lower than that of the forward shock as F^ ax j F max j < 2 x 
10~ 4 K B 1/2 nT , 2 R d . i 7 . Since the typical frequencies v m j and v 1 ^ are around the X-ray 
band in the early afterglow phase, the forward shock emission should peak at a low cooling 
frequency v c j <C v x (this could be achieved with a moderate density ambient medium), 
otherwise the SSC emission is masked by the forward shock emission. 

6. Light Curves of SSC emission 

After a reverse shock crosses the shell, the hydrodynamic quantities of a fluid element in 
the shell evolves as the Lorentz factor T oc R~ 7 / 2 , the mass density p oc R~ 3 T and the internal 
energy e oc p 4 / 3 for a thick shell case, while for a thin shell case they can be approximated 
as T oc R- 2 , p oc (R- 3 T) 6 / 7 and e oc p 4/3 (Kobayashi & Sari 2000). Using these scalings, 
one obtains that the flux below the peak frequency min[z/f , z/ c ] evolves as ~ t" 1 / 2 for a 
thick shell and ~ t~ 2/3 for a thin shell. The flux between and v[ c falls as ~ £-(5p+i)/5 
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for a thick shell and ~ £~( 3 p +1 )/ 3 for a thin shell. When uj, c ~ t~ 2 becomes lower than the 
observed frequency, the flux drops exponentially with time. However, the angular time delay 
effect prevents an abrupt disappearance. The flux will be determined by off-axis emission 
if the line-of-sight emission decays faster than t~ 2+lS where (3 = — (p — l)/2 or — p/2 is the 
spectral index (Kumar & Panaitescu 2000). 

The temporal index of a SSC flare in the rising phase is rather uncertain, because it 
depends on the density profile of the fireball shell. We assume a simple homogeneous shell. 
The hydrodynamic quantities and the number of shocked electrons during shock crossing are 
given by T oc i?" 1 / 2 , e oc T 2 , p oc R' 2 Y^ and N e oc R 2 for a thick shell, and T oc R~ 9 , e oc T 2 
and p oc i? _3 r _1 and N e oc R?l 2 for a thin shell where g is a parameter (e.g. Kobayashi 
2000). Around the shock crossing, the shock becomes mildly relativistic (Sari & Piran f995; 
Kobayashi, Piran & Sari 1999). We consider the behavior of a light curve near the peak 
(g = 1/2). The break frequencies of the SSC emission evolve as v* c oc t~ 3 / 2 , i/f oc t 1 / 2 
(thick shell) or t 5 / 2 (thin shell). Since during the shock crossing z/f increases while v]p 
decreases, an interesting regime is v 1 ^ < u b s < v]F . Assuming p = 2.5, for a thick shell the 
flux increases as t 7 ^ 8 in this regime, and it is almost constant in the other regimes. For a thin 
shell, the flux in this regime increases as t 15//8 . If z/^p (or v^p) crosses u b s , the flux behaves 
as t~ 5 / 6 (or t 9 / 8 ). If both break frequencies pass the observed frequency v]p < v obs < , 
the flux decays as t™ 3 / 4 . 



7. X-ray Flare from Thin Shell 

In this section, we first study the typical thin shell case, and show the spectrum and light 
curve of the early afterglow. If the reverse shock is in the thin shell regime, the SSC peak is 
separated from the GRB emission. The deceleration time is t d ~ 190 C^52 /3 (ro/80)- 8 / 3 (n/5)- 1 / 3 
sec. The typical and cooling frequencies of the forward shock synchrotron emission at the 
deceleration time are given (e.g. Sari, Piran & Narayan 1998) by 
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where we assumed a larger value of e# for the fireball than that of the blast wave. Using eq. 
(4), the break frequencies of the reverse shock SSC emission are 

We plot the broad band spectrum in figure 1. The reverse shock SSC emission (thick solid 
line) dominates X-ray band (5keV ~ 10 18 Hz) in this example case. Figure 2 shows an 
X-ray light curve around the deceleration time. The reverse shock SSC emission produces 
an X-ray flare around the deceleration time. For reference, we plot the X-ray light curve 
of GRB 050406. Compared to the the theoretical SSC flare (thick solid line), the observed 
flare rises more rapidly, especially around the peak. As we gave a caveat in the previous 
section, the temporal index in the rising phase is rather uncertain. This sharp rise might be 
due to inhomogeneity of a fireball shell (we have assumed a homogeneous shell to evaluate 
the light curve). The internal shock model requires a highly irregular outflow from the GRB 
central engine. Since the hydrodynamic interaction inside the flow smooths the velocity and 
pressure profiles, but not the density profile, fireball ejecta might have an irregular density 
profile at the deceleration time. Emission from the ejecta during a reverse shock crossing 
could reflect the light curve of the prompt emission produced by internal shocks (Nakar & 
Piran 2004). The duration of GRB 050406 was T 90 = 5 ± 1 s in the 15-350keV band and 
the light curve peak had a fast rise, exponential decay (FRED) profile (Krimm et al. 2005). 
Around the peak, a reverse shock might hit a higher density part of the shell. 

Another possibility is that a reverse shock might stay in the Newtonian regime during 
the whole evolution. The SSC emission from a Newtonian reverse shock is expected to rise 
as rapidly as oc t 4p ~ 2 = t 8 for p = 2.5 (thick dashed line). In such a case, the forward shock 
emission also increases faster (thin dashed line). This could happen if later ejecta from the 
central engine injects additional energy to the inner tail of the shell after the internal shock 
phase, and it further smooths the velocity profile of the shell. When the spreading of a 
shell width is not significant, a reverse shock does not evolve to mildly relativistic at the 
deceleration time (Sari & Piran 1995). X-ray flares with a moderate amplitude like this event 
can be produced by the reverse shock SSC process, although the sharp structure dt/t < 1 
might be difficult to be explained in external shock related models (Kobayashi & Zhang 
2006). 

We consider how the relative amplitude of X-ray flares depends on parameters. Since 
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X-ray flares occur a few hundred seconds after the prompt emission, we will consider the thin 
shell case. Fixing the deceleration time t d ~ 190 C^5 1 2 /3 (r /80)^ 8 / 3 (n/5)- 1 / 3 sec, the initial 
Lorentz factor is a function of the ISM density (we assume the typical explosion energy and 
redshift). A larger ratio of e e /e# enhances the scattering process. Assuming the equipartition 
value e e = 0.3, we evaluate the relative amplitude of an X-ray flare as a function of e#. The 
results are shown for different values of n in figure 3. The amplitude is as large as a few tens 
if e B <C 1 and n>l. However, a high density requires a low Lorentz factor (i.e. T ~ 55 for 
n = 100). 



8. Conclusions 

We have investigated the synchrotron self inverse-Compton (SSC) emission from the 
reverse shock. The synchrotron process is expected to produce optical/IR photons, which 
are up-scattered into the X-ray band by electrons heated by the reverse shock. For a thick 
shell, the X-ray flare occurs at the end of the prompt gamma-ray phase. The emission decays 
as ~ t~ 2 or slightly steeper. For a thin shell, the emission initially increases as ~ t 2 (the 
scaling could be significantly different if a fireball shell is highly irregular), the peak should 
be separated from the prompt emission, and after the peak the flux decays as ~ £~( 3 p +1 )/ 3 or 
steeper. If off- axis radiation dominates, the temporal and spectral indices L oc t'V' 3 should 
satisfy the relation a = —2 + f3. A weakly magnetized fireball (e#,/ < e# <C 1) in a high 
density ambient medium provides favorable conditions for producing a significant X-ray flare. 
The contrast between the IC flare and the baseline X-ray emission is at most one order of 
magnitude if the synchrotron process dominates the electron cooling. If the IC dominates, 
the contrast could be larger. However, since most of the energy is radiated in the 2nd order 
IC component around 10-100 MeV, the X-ray bump and the baseline X-ray emission are less 
energetic than the latter. The optical flash (due to synchrotron) is highly suppressed. 

Recently Swift XRT detected X-ray flares in the early afterglows of GRB 050406 and 
GRB 050502b (Burrows et al 2005). The afterglow of GRB 050406 brightens by a factor of 6 
between 100 and 200 s post-burst before starting on the rapid decline seen in other prompt 
X-ray afterglows. GRB 050502b had an even stronger X-ray flare, brightening by a factor of 
~ 10 3 to a peak 700 s after the burst. Both afterglows were very faint at 100 s post-burst (a 
factor of a few - 100 fainter than previous XRT-detected afterglows). 

In the case of GRB 050502b, the X-ray flare contrast factor ~ 10 3 requires a Compton 
parameter x larger than ~ 100, which means e e /es > 10 8 . With plausible values of the 
other parameters, the SSC emission should appear at an energy band well above X-rays. 
The very sharp profile dt/t <C 1 also rules out external shock models. Another mechanism 
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(e.g late time internal shocks; Burrows et al. 2005; Falcone et al. 2006; Romano et al 2006; 
Zhang et al. 2006; Nousek et al. 2006; Fan & Wei 2005; Wu et al. 2005; Ioka, Kobayashi & 
Zhang 2005) is likely to play a role in the production of the X-ray flare. On the other hand, 
the X-ray flare of GRB 050406 might be explained with the reverse shock SSC emission 
discussed here, although the rather narrow feature dt/t < 1 might disfavor the external 
shock interpretation. If we apply the SSC model to the flare, the peak time td ~ 200 sec 
gives an initial Lorentz factor r ~ 80 (t (i /200s) 3 / 8 C 3/8 J B 5 1 2 /8 (n/5)- 1 / 8 . The reverse shock SSC 
emission can not explain multiple X-ray flares in a single event. Such behavior is observed in 
recent Swift bursts (Burrows et al. 2005; Falcone et al. 2006; Romano et al. 2006; O'Brien 
et al. 2006). 

In thick shell cases, the SSC flare rising portion overlaps in time with the prompt 
emission seen by BAT. Although it could be difficult to separate the X-ray emission of the 
prompt and the reverse SSC flare rising components, the rapidly decaying portion of the 
reverse SSC flare could be detectable. Swift has reported many of the X-ray afterglows 
detected have an early steep decay phase, before going into a more common shallow decay. 
In some of these cases, the early steep decay may be interpreted as the tail portion of a 
reverse shock SSC X-ray flare. 

We note also that the reverse shock SSC mechanism predicts, besides an X-ray flare, 
a strong GeV flare from 2nd order IC, and in some cases from first order IC. Thus, bursts 
with strong early X-ray flares should be good candidates for GLAST. 

We thank Patrizia Romano and the Swift team for providing XRT data. This work 
is supported by Eberly Research Funds of Penn State and by the Center for Gravitational 
Wave Physics funded by NSF under cooperative agreement PHY 01-14375 (for SK), NASA 
NNG04GD51G (for BZ), NASA AST 0098416 and NASA NAG5-13286 (for PM), and NASA 
Swift GI program (for BZ, SK and PM). 
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Fig. I. — Spectrum at the deceleration time: Reverse Shock SSC (thick solid) and Syn- 
chrotron (thin dashed) emission, and Forward Shock Synchrotron emission (thin solid). 
z = l,E 52 = l,r = 80, n = 5,e e = 0.3, e B = 0.03 and e B j = 0.01 The flux is normal- 
ized at the forward shock peak. 
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Fig. 2. — X-ray Light Curve: Reverse Shock SSC (thick solid) and Forward Shock Syn- 
chrotron emission (thin solid). The parameters are the same as in figure 1 except Tq = 79. 
The circles represent the X-ray light curve (counts s _1 ) of GRB 050406 from Romano et al. 
(2006). The theoretical light curves are normalized as the reverse shock SSC emission at the 
deceleration time fits the observed peak (the peak counts ~ 6 counts s _1 ). If the width of 
the fireball shell is constant before the deceleration, the shock emissions rise more rapidly at 
the beginning (dashed lines). 
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Fig. 3. — X-ray flare amplitude: Flux ratio between the reverse shock IC and forward shock 
synchrotron emission Fl° / F v j at X-ray band (5keV) at the deceleration time is plotted as a 
function of 6b- n = 5 (solid), 10 (dashed), 50 (dashed dotted) or 100 (dotted), z = l,E 52 = 
l,e e = 0.3 and £b/£bj — 3. For a smaller 6b, the Compton parameter x is larger and the 
typical frequency of the IC emission i/p is lower. The amplitude of an X-ray flare peaks at 
a moderate e# with which i/p is close to the X-ray band. 



